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Abstract 

J>->, We propose a new method to study motions of mixtures in fluid interfaces. We 

extend the equations of equihbrium in interfaces and the results associated with 

^ ■ travehng waves for van der Waals hke fluids [21]. Maxwell rule is extended to in- 

^+H , terfaces of fluid mixtures out of equilibrium. Formulae like Clapeyron relation are 

Y^ \ obtained for isothermal layers. 
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00 . 1 Introduction 

o 
> 

KA , In this paper we propose a new attempt to study motions in interfaces of fluid 

H \ mixtures. First, we recall the study of interfaces between mixture bulks in 

equilibrium. We remind problems for an interface moving in a single fluid and 
extend them to cases of mixtures. 

How to calculate the proportions of constituents in each bulk of a n-components 
mixture? This problem is solved in the case of isothermal equilibrium for mix- 
tures whose free energy is known [7]. A simple calculus consists in writing the 
total free energy of the mixture is minimum. 

Practically, calculus are available only if: 



Email address: henri.gouin@univ-cezanne.fr (Henri Gouin). 
Preprint submitted to : IMA Volumes in Mathematics and its Applications, 52, pp. 111-121 (1993). 



- the container has a simple shape allowing one dimension calculus (cylinder 
or sphere) 

- the mixture has a few number of components 

- the mixture free energy is not a too complex function of component densities 
and temperature. 

n 

Let us denote by p^ the density of constituent «; p = X! Pi '^^ ^^^ density 

of the mixture, v = — the specific volume and Cj = — the concentration 

P P 

of constituent i. The volume free energy at temperature 9 is in the form 

Go = Go{v,Ci,9). The total free energy is the sum of the thermodynamic 

free energy and the interfacial energy. 

The determination of the function Go is certainly of first importance, but it 
is not the topic of this paper. This is both a theoretical and an experimental 
problem and an equation of state with coefficients given by experiments is 
used for calculations. The free energy must be independent of the proportion 
of constituents when the density is vanishing [5,15,18]. 

In the case of flat interface, the knowledge of interfacial energy is not neces- 
sary for the calculation of balance equations between the bulks. The knowl- 
edge of interfacial energy is necessary to study spherical interfaces (bulks, 
drops, aerosols...). A Landau- Ginzburg model consisting in a quadratic form 
of density gradients representing the interfacial energy [12]. Coefficients of the 
quadratic form are assumed to be constant (in the case of a single fluid, this 
corresponds to a mean field theory [16]). A quadratic form of the gradients of 
speciflc volumes can be also used; nevertheless, assuming that the coefficients 
of the quadratic form are constant, it has a different physical meaning. In fact, 
the choice of a model is of no consequence outside of the interfacial layer. The 
Landau- Ginzburg model represents an interface and the bulks and builds a 
complete theory of the mixtures in dynamics [8,9,13]. 

In this paper, we apply the model to motions through interfaces for both single 

fluids and mixtures. Because of physic scales, motions can be often considered 

as permanent [4]. Interface moving in a single fluid was investigated by Slemrod 

[20-21] as a traveling wave whose image in a reference space is propagating with 

the constant velocity C. It is equivalent to say the fluid crosses a stationary 

interface with the flow q [4] (g is the product of velocity and density). When 

dissipative phenomena are neglected, the problem boils down to equilibrium 

X 
with a speciflc free energy increased with — q^ v^. 

The general study of mixtures in dynamics can be applied in the case of 



interfaces [9]. The constituent i crosses the interface with the flow g^. We are 
back to a static problem with the same specific free energy and an additional 
term in the form 

Each constituent of the mixture has its own reference space [1,9,17]. Due to 
the fact that the interface has different velocities in each reference space, the 
generalization of traveling wave is less intuitive. 

For a n-component mixture, a permanent motion through an interface yields 
n + 1 equations of dynamical equilibrium. The equations represent a linear 
system of the n + 1 variables q]^ and 9. When the specific volume v and the 
concentrations Cj are given in each bulk, a simple calculus yields the temper- 
ature of the mixture and the flow for each constituent. Obviously, solutions 
are acceptable only if values q^ and 9 are positive. 

Surprisingly, the problem of an interface in equilibrium is more difficult to 
study than the dynamical problem. It is more difficult to choose g^^ and 9 {qi 
is null in equilibrium case) and to deduce concentrations and specific volume 
than to choose concentrations and specific volume in each bulk and to deduce 
qf and 9 from a system of linear equations. 



2 Isothermal motion of a single fluid in the "Korteweg - van der 
Waals" theory of capillarity. 



Let us consider the one dimensional motion of van der Waals like fiuid. We 
notice that the same equations describe a traveling wave of dynamic phase 
transition and a fiow through an interface in the case of permanent motion. 

M. Slemrod [21] consider the one-dimensional motion of fiuid possessing a 
specific free energy of the forno: 

G{v,9)=G^{v,9) + ^{^r 

where X is the lagrangian mass variable, e is a small positive parameter. In 
the referential space, the balance laws of motion may be written in the form: 



^ The notations are different from the ones given in [20-21] 
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where V denotes the velocity of the fluid, /x the viscosity, p the pressure, E 
the specific total energy, k the coefficient proportionality between the heat 
flux and the gradient of temperature in the Fourier Law and ^ = X — Ct. 
The traveling wave theory of phase transitions attempts to determine when 
two homogeneous phases may be joined by a traveling wave solution {v = 
v{X - Ct), V = vlx - Ct), q = q{X- Ct)} of equations (1). Slemrod flnds 
that the volume mass is solution of the equation: 



'— = - C^(^y - v^) - (^p - p^) - ^c V 



(2) 



Po and Vo denotes the values of p and v in the bulks. Because of the physical 
scales far from critical conditions, we established in [4] the motions of van der 
Waals fluids crossing interlayers can be considered as permanent. In the case 
of permanent motion of inviscid fluid, Eq. (1)2 yields 



dV 
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d V 
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where x denotes eulerian variable and 7- — = v 20 . From V = qv where q 

oX 
denotes the constant flow of the fluid for an unidimensional permanent motion, 

we deduce: 



ax2 



-q {v - Vo) - (p - Po) 



At t given, the form of the equation is similar to Eq. (2), which allows to 
obtain interpretations about dynamical Maxwell rule as in [20-21]. 



3 Isothermal motion of fluid mixture in the "Korteweg 
Waals" theory of capillarity [9] 



van der 



For the sake of simplicity, we study a mixture of two fluids. The method can 
be immediately extended to any number of constituents. No assumption has 
to be done about composition or miscibility. 



The motion of a two- fluid continuum can be represented with two surjective 
differentiable mappings: 

z^ Xi = Mi{z) and ^ ^ X2 = M2{z) 

(Subscripts 1 and 2 are associated with each constituent of the mixture) 

Xi and X2 denote the positions of each constituent in reference spaces Dqi 
and Dq2- The lagrangian of the mixture is: 

1 1 

^=2 ^^^^ ^ 2 ^^^^ ~^~ ^^^^ ~ ^^^^ 

where Vi and V2 denote the velocity vectors of each constituent, p^ and pg 
are the densities, f2i and ^^2 are the extraneous force potentials depending 
only on z = (t, x) and e is the volume internal energy. 

The expansion of the Lagrangian is in a general form. In fact dissipative phe- 
nomena imply that Vi is almost equal to V2 and we do not take into account 
some kinetic energy associated with the relative velocity of the components. 
Because of the interaction between the constituents, e does not divide into 
energies related to each constituent of the mixture, like for simple mixtures of 
Huids [14]. The mixture is supposed to be no chemically reacting. Conserva- 
tions of masses require: 

p,detF, = p„,TO (3) 

where subscript i belongs to {1,2}. At t fixed, the Jacobian associated with 
Mi is denoted by Fj and p^j is the reference specific mass in Dgi- 

In differentiable cases, Eq. (3) is equivalent to: 

^ + divp,V, = 

The volumic internal energy e is given by the thermodynamic behavior of the 
mixture. Each constituent has a specific mass; in the same way, two specific 
entropies Si and S2 are supposed to be associated with constituents 1 and 2. 

For an internal energy depending on gradients of densities, the volume internal 
energy is: 

e = e{si, S2, Pi, p2, gradpi, gradpa) (4) 

The quantity 

defines the specific enthalpy of the constituent i of the mixture. Subscript 
7 corresponds to the spatial derivatives associated with gradient terms; as 



usually, summation is made on repeated subscript 7. 



9, 



Pi 



deGnes the temperature of the constituent i of the mixture. 
In practise we use the expression 



e = a{si,S2,pi,p2) + 2 ^ 



where Q is a quadratic form with constant coefficients: 

Q = C\ (gradpi)^ + 2L'gradpi gradp2 + C2 (gradp2)^ 

and a{si, S2, Pi, P2) is the value of internal energy in the homogeneous bulks. 

To obtain the equations of motions, we used variational principle whose orig- 
inal feature is to choice variations in reference spaces (Fig. 1) [9]. 




Figure 1 

Variations 61 and 62 of motions of particles are deduced from Xi = iJi^x, t, P^) 
and X2 = 'ip2{x^t, 132)- They are associated with a two-parameter family of 
virtual motions of the mixture. The real motion corresponds to /3^ = and 

Obviously, a classical variation of function in space W (as Serrin's p. 145, [19]) 
with only one parameter cannot give as many informations as the ones given 
by any two-parameters family of virtual displacements. 



The equation of the motion of each constituent of the mixture is: 

Ti = 9i grad Si - grad {hi + Qi) (5) 

Let us note that the motion of each constituent is described by the same equa- 
tion that for a single fluid [10,19]. In applications, the motions are supposed 
to be isothermal and O-i = {9 common temperature value for all the compo- 
nents). This case corresponds to strong heat exchange between components. 

Eq. (5) yields: 

Ti + grad (v^j + VLi) = Q 

where (Pi = hi — 9si is the free specific enthalpy (or chemical potential) of 
constituent i. In the bulks, the volumic free energy associated with a is denoted 
by go = doiPi, P21 ^)- A complete study of the barycentric motion of the mixture 
and of the equation of energy is given in [9] . This generalizes results for single 
fluid [3,6,11]. 



4 Motion of an isothermal fluid mixture through a plane interface 

An interface in a two-phase mixture is generally schematized by a surface 
without thickness. Far from critical conditions, this layer is of molecular size 
and density and entropy gradients are very large. A continuous model schema- 
tizes such areas by using an energy in form (4) extending forms given in [2,16] 
for compressible fluids. 

For the same reasons than for a single fluid, motion is supposed to be isother- 
mal, stationary and one dimensional with respect to x. 

The acceleration of component i is: 

1 d o^ 
' 2dx^pr' 

and the equations of motion yield: 

Fi = grad {Ci Api + D^p^-g . } 



F2 = grad {D Api + C2 Ap^ - g „, } 
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or: 



A Pi 



Dp'; + c2P^' = <?o„, + ^4 + ^2 

^ P2 



(The derivatives along the motion axis are denoted by '). 
Combination of the last two equations yields the first integral: 

7 7 -'- ?i ^ ^.1 / ^ ^~i '2 T-. ' ' -'- ^-f '2\ 7 

go + hpi + k2p2- ^ 2 ''2 ^ ^^ +Dpip2 + ^C2pi) = k3 

In each bulk, densities have zero-gradients. Eliminating constants ki, dynam- 
ical conditions through the interfacial layer yield: 



kpi(pl>P2)] =0 
kp2(Pl'P2)] =0 

[9- Pi9,p^ - P29.P2] 



(6) 



where [ ] denotes the discontinuities through the layer and 

1 ql 1 ql 
9 = 9o~- -— • 

2 Pi 2 P2 

In case of equilibrium {qi and g2 are null), the minimum of the total free 

energy, with a given total mass for each constituent, yields conditions (6). 

1 ^2 ^ ^2 

By adding the term to go, the study of flat interfaces crossed 

2 Pi 2 P2 

by components of a mixture turns back to an equilibrium problem. In fact, a 

complete study of the thickness of the interfacial layer or of spherical interfaces 
of microscopic size requires the non-linear model. 

It is more classical to use the variables v or c. The mapping (p^, p2, g) 



(f , c, G) where v 
allows to write 



1 



P1+P2 



P2 
P 



) Go 



9o 

p 



and G 



Go--v^ 



-c c' 



[G,„(t;,c)] = 

[G,e(t;,c)] = (7) 

[G-vG,^-cG^c] =0 

Conditions (6) or (7) express that the points corresponding to the bulks of 
the Gibbs surface (S) associated with the dynamical free energy g 01 G are 



contact points of a bitangent plane (Fig. 2) 



o + (a) 




Figure 2 

The pressure of the mixture is p = —Go,v{v,c), the dynamical pressure is 
P = -G,,{v,c). 

The ^-coordinate of the point of intersection of the bitangent plane with the 
g-axis is —P . The chemical potential G^c is denoted by $. 

Contact points of the bitangent plane with surface (E) generate a curve. Along 
this curve, the straight line (D) connecting the contact points is the charac- 
teristic line of the plane. Along the line, the relation 

{v — vq) d G^y + (c — Co) d G^c = 

allows to generalize Clapeyron formula such that: 

dP ci — Co 
(i$ f 1 — f 

where subscripts and 1 are associated with the two bulks. 



5 Maxwell rules for fluid mixtures 



At a given temperature 6, the differential form dG = — P dv + ^ dc yields the 
relation: 

/ ($ - $0) dc-(P- Po) dv = 

J{C) 

where (C) is an arbitrary curve connecting the two bulks in the space (c, -u). 
The common values of $ and P in the bulks are denoted by <l>o and Po- 

We can calculate the integral along special paths: 



Along the curve (C) given by the imphcit function c{v) defined such that 
G,c (c(f ), v) = $0 5 the relation 



f (P- Po) dv = 

J{C) 



represents the Mawxell equal area rule for the path abed in the plan (P, v) 
wher^fj P = —G^v (c (v) , v) 

Likely, for P = Pq, the relation: 



/ ($ - $0) dc = 

J{C) 



represents the Mawxell equal area rule for the path abed in the plan (c, $) 
where $ = G^c {c,v{c)) is such that v is expressed as a function of c by the 
relation G,„ {c,v{c)) = — Pq- 




Figure 3 

These results generalize the ones obtained by Slenirod in the case of inviscid 
fiows of single fiuids [21]. The knowledge of the free specific energy Go is 
sufficient to solve the previous problem. Function Go is obtained by integration 
of the equation of state p = p{v, c, 9): 



dGoiv,c,e) 
dv 



-P (v, c, 6) 



Function p is determined except for an additional function of c and 6. With 
the physical assumption that Go is independent of c when v vanishes, the 
additional function depends only on 6. 



(C) is not the curve associated with P = P{c, v) at c constant. 
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The simplest model of equation of state is by van der Waals: 



R9 _ a 

V — b v"^ 



P=^^.--, (9) 



where a and h are given by the mixing rules: 

a^ = {I- c)al + cal et 6 = (1 - c)bi + 062 (10) 

(Subscripts 1 and 2 are associated with the components 1 and 2). 

One deduces: Go = —R9Ln (v — b) . 

V 

Other equations of state are proposed for example p = — (par- 
ti — 6 v{v + b) 

ticularly in petroleum industry [5,15,18]) or by changing the form of the mixing 

n n 

rules (other expressions than a^ = "^ of Ci and 6 = ^ 6i Cj can be consid- 
ered) . 

Equations (7), (8), (9) and (10) lead to the system: 

e V 2 V 2 (b+l)\ 



V — b 1 — c c v'^ 

9{l - Ln{v -b) -\ -\ + - — —Qi - 2- 



v-b' 2 (l-cy ' V 

_^l,i!L^2 ,(&2 + l)(&+l) 

v-b^ 2 c2 






which is a linear system with respect to 6, qf, g2- 



6 Conclusion 



The study of fluid mixture motions by a model taking into account density 
gradients of components requires the knowledge of the global free energy. In 
fact, only physical experiments and molecular theories provide an equation of 
state (in van der Waals type). Isothermal motions can be studied by using 
such an equation and yield the jump conditions through interfacial layers. To 
carry the program one step further and to study non-isothermal motions, it is 
necessary to get additional knowledges (such that specific heats...). 
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